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ABSTRACT 

An analytical model to study the coupled transverse and longitudinal vibrations of a single lap adhesive joint is 
proposed in this paper which includes partial differential form of the motion equations. A balanced single lap adhesive 
joint consist of two identical adherents of mild steel which are lap jointed over a certain length by a viscoelastic material, 
epoxy resin(araldite). Adherents are modeled as Euler-Bernoulli Free-Free beam. Both transverse and axial deformation of 
adherents, shear and peel stresses at the adhesive joint interface and deflection of mid plane of adhesive layer (3 -parameter 
model) is considered in deriving the equations of motion. The classical two parameter elastic foundation model violates the 
equilibrium condition of the adhesive layer; to eliminate this flaw, a new three parameter elastic foundation model is 
considered which satisfies the equilibrium condition of the adhesive layer. The governing equations of motions are derived 
for three parameter elastic foundation model. The numerical solutions of the governing equations for free vibrations yield 
the system natural frequency. Experimentation carried out on both monolithic and adhesively jointed beam to observe the 
effect of joint. The effect of thickness of adhesive and joint overlap ratio on system natural frequencies are investigated. 

KEYWORDS: Adhesive Joint, Euler-Bernoulli Beam, Natural Frequency, Monolithic Beam 

INTRODUCTION 

Joining components by using adhesives is becoming more popular with the development of adhesives with high 
adhesion properties. In many structures such as those for flight, automobiles and space vehicles etc. adhesively bonded 
structures have often been used recently because of the great advances in adhesive bonding techniques. Single lap joints are 
used extensively for adhesively bonded test specimens and in various configurations for wide range of industrial products. 

The paper by Goland and Reissner [1] is regarded as a classical work in the area of static analysis of a single lap 
adhesive joint. They derived equations for evaluating the shearing and normal stresses in the bond layer; peel and shear 
stresses were assumed to be constants across the adhesive thickness. A comprehensive literature survey on this subject has 
been presented in [2, 3, 4, 5, 12, 13, 15]. Saito and Tani [6] have derived equation for predicting the modal parameters of 
the coupled longitudinal and flexural vibrations of a system consisting of a pair of elastic beams fixed at both ends; the 
longitudinal deformation of the adhesive layer is neglected in deriving the equations of motion. He and Rao [8, 9] 
developed an analytical model to study vibrations of bonded joints subjected to transverse loadings. The unjointed ends of 
the beams were assumed to be simply supported. The equation of motion was developed using Hamilton's principle. Vaziri 
A, N. Hashemi [10] developed a model to understand the effects of defects such as a void in the overlap on the system 
response. As described in [6, 8, 9, 10], the peel stresses acting on top and bottom surfaces of adhesive layer are assumed 
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identical, which violates the equilibrium condition of the adhesive layer. The paper [16] deals with free vibration analyses 
of adhesively bonded plates. The eight-node isoparametric plate element was used to model the adherends, they show the 
influence of patch size on the natural frequencies and mode shapes of the five -layer cantilever bonded plate with a stepped 
type. A three-dimensional free vibration and stress analysis of an adhesively bonded functionally graded single lap joint 
was carried out by [17]. They show that the effect of the adhesive material properties, were negligible on the first ten 
natural frequencies and mode shapes of the adhesive joint. Both the finite element method and the back propagation 
artificial neural network (ANN) method were used to investigate the effects of the geometrical parameters. The influence 
of adhesive properties on the transverse free vibration characteristics of single lap jointed cantilevered beams were studied 
using FEM [18], they show that the transverse natural frequencies of beam increase with increase in Young's modulus of 
adhesive where as poisson's ratio has little effect on frequencies. T-C Ko, C-C Lin, R-C Chu [19] used an isoparametric 
adhesive interface element in the vibration analysis of adhesively bonded structures. In their investigation, the natural 
frequencies of laminated plates with a lap joint were found to be reasonable when compared with the results for plates 
without joints. Gang Li, P. Lee-Sullivan [20] used a geometrically nonlinear two-dimensional FE method to study the 
adhesively bonded balanced single-lap joints. The effects of plane strain and plane stress conditions, boundary conditions 
at the adherent ends, filleted and unfilleted overlap end geometries and two different adhesive materials on the bending 
moment factor and the adhesive stresses have been evaluated. A three-parameter, elastic foundation model is proposed by 
[21, 22] to analyze interface stresses of adhesively bonded joints 

In this study a novel three parameter elastic foundation model of adhesive joint is considered. The objective of the 
paper is to examine the effects of single lap adhesive joint on the modal parameters of the structure. A deflection of mid 
plane of adhesive layer (3 -parameter model) is considered in deriving the equations of motion; the classical two parameter 
elastic foundation model violates the equilibrium condition of the adhesive layer; to eliminate this flaw, a new three 
parameter elastic foundation model is considered which satisfies the equilibrium condition of the adhesive layer. A beam 
with Free-Free configuration is considered in analytical and experimental investigation. 

THEORETICAL INVESTIGATION 
Single Lap Adhesive Joint 

For the vibration analysis of adhesively bonded lap joint, the following assumptions are made: (1) The beams are 
purely elastic and the adhesive is linearly viscoelastic. (2) The system deforms in plain stress or plain strain and the usual 
Euler-Bernoulli beam theory is considered. (3) Shear stress in the bond material varies linearly through the adhesive 
thickness. (4) The adhesive layer usually has much lower extensional and bending stiffness than those of the adherents; 
therefore the axial and bending moment of the adhesive layer can be ignored. 

Schematic diagram of a balanced single lap adhesive joint is shown in Figure 1. The mild steel adherents are 
modeled as Euler Bernoulli beams supported on viscoelastic foundation, which resist both peeling and shear deformations. 
The beams are purely elastic and have thickness ti, width w, length (1+c), density p and Young's modulus E. The adhesive 
is considered to be viscoelastic and its complex shear and elastic modulus is assumed to be G a and E a respectively, h 0 is the 
adhesive thickness. 
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Adherent 1,(E,I,A) 



Region 3 



Adherent 2,(E,I,A) 



Figure 1: Schematic Model for a Single Lap Adhesive Joint 
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Figure 2: Free Body Diagram of an Element in Region 2 

Considering the free body diagram of an element in region 2 (jointed region), (Figure 2), the equations of motion 
for each beam element are: 

For Adherent 1 



dx 



d 2 y, 
dt 2 



(l) 



_ dM } u 

Q x - + — zw = 0 

dx 2 

T\ EA * \-rw-(pA)— 1 - = 0 

dx\y dx J dt 



(2) 
(3) 



For Adherent 2 



(4) 



<9m 



(5) 



^ f ^ A du 1 \ , A s d 2 u 1 

dx\ dx J dt 2 



(6) 



Where, <2^ is the shear force and M q (q =1,2) the bending moment at the edge of the joint element, A is the cross 
sectional area of the beam, u h u 2 , yi andy 2 are the longitudinal and transverse displacements of adherents. Figure3 shows a 
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three-parameter elastic foundation model for adhesive layer. 
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Figure 3: Three-Parameter Elastic Foundation Model for Adhesive Layer 

In a two parameter model of adhesive layer ; the peel stresses acting on the top and bottom surface of the adhesive 

dz 

layer are assumed identical i.e. o x = o 2 then the equilibrium condition of the adhesive layer gives — = 0 ; this equation 

dx 

suggests that the shear stress is constant along the x-direction which contradicts the shear stress distribution predicted by 
two-parameter elastic foundation model [21]. To remove this contradiction we can assume that <Jx ^ °" 2 
Then we have, 
dx o 2 —G\ 



dx h 0 



*0 



(7) 



Shear stress in the bonded area considering longitudinal and transverse displacement of adherent and deflection y a 
of mid plane of adhesive layer can be expresses as [21]. 



2/Zn 



2u x - 2u 2 + t x 



dx dx 



+ G, 



dx 



(8) 



And the adhesive peel stresses can be expresses as [21] 



vi=—r—(yi-ya) ; °i =^^(y a -yi) 



(9) 



The complex shear modulus G a = G'(l+u.ag) and complex elastic modulus E a = E'(l+Tj ae ), where G' and E' are the 
storage moduli and r| ag and r| ae are the adhesive loss factors. 

Longitudinal forces and moments in the beams from beam theory are - 



N q =EA(du q /dx) 



and 



M q =EI(d 2 y q /dx 2 ) ,q=l,2 



(10) 



Now, from equations 3 and 8 we have 



EA 



( '^2 ^ 

o u 



dx 2 



-(pA) 
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2K 
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(11) 



Combining equations 1, 2, 9 andlO and simplifying gives, 
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Similarly, from equations 6 and 8 we have 
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Combining equations 4,5, 9 and 10 and simplifying gives, 
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Considering the equilibrium of adhesive layer, equation 7, 8 and 9 gives 
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(15) 



Now assuming that the beam is excited with initial harmonic excitation of frequency co, the displacement field is 
Uj(x,t)=Uj(x). e i(Dt and y r (x,t)=Y r (x). e io3t ; j=l,2 and r=l,2,a. The equations 11 to 15 becomes - 
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(20) 



In the usual fashion above equations (16 to 20) are transformed into five algebraic equations by assuming 
Uj=Z n .e Kx and Y q = z\.e Kx where ; j =1,2 and q=l,2,a (21) 
Where Z n and Z n * are the arbitrary constants to be determined from boundary conditions. The equations (16 to 20) 
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and (21) can be written in matrix form as- 
[J] 5 x 5 {V}=0 



(22) 



where {Y}=[Z l Z 2 Z* Z 2 * Z a *] T and [J] 5 x5=[aij] 5 x5 where, 
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For non-trivial solutions of equation (22) , 
det [J] 5x5 =0 



(23) 



Expansion of equation (23) gives fourteenth degree polynomial in Z, coefficient of which includes known 
geometrical and material properties of adherents and adhesive with only unknown parameter co. The 14 roots of the 
polynomial are A, n ,(n= 1 ,2 .... 1 4) . 



14 14 

Now, Uj =^Z n .e Kx ; Y q =^Z* n .e Kx j =1,2 and q=l,2,a 



(24) 



For unbonded part of the upper and lower beam i.e. region 1 and region 3, the equilibrium equations are those of 
classical Euler-Bernoulli beam equations (Figure 4). The equations of motion of simple beams (unbonded part) for 
transverse and longitudinal vibrations are - 
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Figure 4: Bonded Joint System and Co-Ordinates 
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Assuming that yu(x n ,t) =Y n (x n )e ltot and u n (x n ,t)=U n (x n )e mt .The solution to above equations is in the form of : 

4 4 2 2 

Y n =Y,R n e SnXu ; Y 22 ^ P n eSnX22 ; U n =^J n e D » x " ; U 22 = ^n^ 11 (27) 
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Now with boundary conditions and continuity equations between bonded and unbonded regions of the beam, we 
can obtain frequency equation. 

BOUNDARY CONDITIONS AND FREQUENCY EQUATION FOR SINGLE LAP JOINT 

The boundary conditions for adhesively jointed free-free beam shown in Figure 4, are as follows: 
Boundary conditions for free ends of the beam at x n =0 and x 22 = 1 , 



For adherent 1, at x n =0 : Un=0 ; 
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0 ; 
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0 (Region 1) 



For adherent 2, at x 2 2= 1 : U 2 2=0 ; 
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(28) 



For the free ends of the joint region (Region 2) at x= 0 and x= C 
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The continuity conditions of the beams 
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at x 22 =0 ; x= C : F 22 = Y 2 ; 
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(30) 



Where Q is the shear force at the jointed region, assuming shear stress x at the free end of joint; the shear force Q 
is expresses as - 



2U l - 2U 2 + t x 



dY dY 2 
— L + — - 
dx dx 



dYa_ 
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(31) 



Using equations (24) and (27); the above 26 homogeneous equations, (28) to (30) may be presented in matrix 



form as- 



[T][Z]=0 



(32) 



In which [T] is a [26x26] coefficient matrix and [Z] is a column vector of 26 constants as, Z n (n=l,2....14); R n , P n 
(n=l . . ..4); J n , H n (n=l,2) , for non-trivial solution of equation (32) the determinant of matrix - 
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[T] =0 (33) 

The equation (33) is the frequency equation, the roots of which yield the complex natural frequencies (cor+ icoi) of 
the joint system. The equation will be in terms of mechanical and geometrical properties of adherent and adhesive. 
The solution of equation (33) was obtained by developing code in MATLAB 7.1. 

EXPERIMENTAL INVESTIGATION: RESULTS AND DISCUSSIONS 

A single lap adhesively bonded joint was prepared using mild steel adherents having thickness 5 mm and width 50 
mm and were joined together using epoxy resin (araldite) adhesive of 0.5 mm thickness. The overlap area was (50x100) 
mm 2 and length of beam L was 600 mm. The beam was suspended at each end with a thin nylon cord to simulate the free- 
free boundary condition in transverse vibration. 

Experimentation was carried out on both monolithic beam and single lap adhesively bonded beam. Figure 5 shows 
the test setup for measuring the dynamic response of single lap adhesive joint and monolithic beam. 




Figure 5: The Test Setup 

A monolithic beam i.e. a beam possessing identical geometrical and material properties to the jointed structure but 
with no joint interface is shown in Figure 6, 



1 


c 


1 * ' w 









Figure 6: Monolithic Beam 



Dynamic response of the bonded joint and monolithic beam was obtained by striking the beam with nylon 
hammer and measuring the response using accelerometer, Figure 7 shows the frequency spectrum of a bonded joint used 
to identify resonant frequencies; commercial package NI-LABVIEW was used to acquire response signals through 
Ni-DAQmx. The natural frequencies thus obtained are compared in Table 3. This preliminary experimentation was carried 
out to show that; adhesive joint in structure does affect the modal parameters of beam. 

Detailed theoretical and experimental investigation was then carried out for single lap adhesive joint. 
The mechanical and geometrical properties of adhesive and adherents used in the analytical and experimental investigation 
are listed in Table 1 and 2 respectively. 
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Table 1: Mechanical Properties for Adherent and Adhesive 
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E' 


G' 


P (Kg/m 3 ) 




Adherent (Mild steel) 


207 GPa 






7850 




Adhesive (Epoxy-araldite) 




2.46 Gpa 


0.83 
GPa 


1150 


0.1 



Table 2: Geometrical Properties of the Model 
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Figure 7: Frequency Spectrum of a Bonded Joint, Used to Identify Resonant Frequencies 



Table 3: Natural Frequencies for Adhesive Joint and Monolithic Beam 



Natural 
Frequencies (Hz) 


Adhesively Jointed Beam 
(Experimental) 


Monolithic Beam 
(Experimental) 


Monolithic Beam 

(FEM) 


COi 


80.795 


84.76 


84.746 


co 2 


201.325 


202.649 


202.85 




431.788 


458.278 


449.23 




660.68 


665.954 


666.22 


co 5 






1138.5 



The natural frequencies co R determined analytically using equation (33) for different overlap ratio (c/L). The total 
length L of the beam was kept constant as the lap ratio is increased. The variation of the normalized resonance frequencies 
<D*(= cdr/cdq) with overlap ratio c/L is shown in Figure 8 ; where co 0 is the fundamental natural frequency of Free-Free beam 
with same dimensions as beam shown in Figure 6. Three more samples of single lap adhesive joint were prepared 
(in addition to one for which geometrical properties are given in Table 2) with c/L=0. 1,0.2 and 0.4 with same thickness of 
adhesive h 0 =0.5 mm. The length L was kept constant equal to 600 mm for all samples. 




Joint Overlap Ratio (c/L) Joint Overlap Ratio (c/L) 

Figure 8: Variation of Normalized Resonance Frequency of Joint System with Overlap Ratio 

It is evident from Figure 8 that there is increase in resonance frequency of adhesive joint system with overlap 
ratio. The results seem to be consistent physically; an increase in lap ratio will make the joint stiffer, thus increase the 
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resonance frequency of the joint system. The theoretical results matches closely with experimental results, it is clear that 
the theoretical model derived here can be used to predict modal parameters of the adhesive joint system. 

Four more samples of beam with single lap adhesive joint were prepared for experimentation having geometrical 
parameters as given in Table 2 and with different adhesive thickness h 0 =0.5, 1.0, 2.0 and 3.0 mm. Natural frequencies are 
measured for each sample and also determined analytically. Figure 9 shows variation of normalized natural frequencies go* 
with adhesive thickness h 0 . There is decrease in normalized resonance frequency with increase in adhesive thickness 
indicating that the joint is becoming more flexible; the joint system stiffness reduces with increasing thickness of adhesive. 
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Figure 9: Effect of Adhesive Thickness on System Resonance Frequency for C/L=0.167 
CONCLUSIONS 

A three -parameter model of adhesive joint was investigated theoretically and experimentally to describe the 
behavior of the joint system. A so called third parameter y a Le . the deflection of mid-plane of adhesive layer was introduced 
in theoretical derivation of the model. Experimentation was carried out on Free-Free Euler-Bernoulli beam with single lap 
adhesive joint having different overlap geometry and adhesive thickness. The results of the experimentation compared with 
those obtained from theoretical model. Effect of joint overlap ratio and thickness of adhesive on resonance frequencies was 
investigated and presented in this paper. 

There is increase in resonance frequency of adhesive joint system with overlap ratio. An increase in lap ratio will 
make the joint stiff er, thus increase the resonance frequency of the joint system. The error, particularly in higher mode is 
more due to systematic error in the measurement of higher natural frequencies. 

The effect of adhesive thickness on the system natural frequencies presented here indicates that the system natural 
frequencies decrease with increasing adhesive thickness. This can be explained considering the system stiffness is reduced 
with increasing the adhesive thickness. 
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